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A. Purpose/Description: This program solves sets of non-homogeneous 
simultaneous linear equations and provides either printed or punched output. 
It is designed for ease of use. Operating instructions and error messages 
are automatically typed. Data values are entered in free form notation as a 
group of digits with a decimal point. An optional power of ten may be added 
to each value. 

B. Method : The Jordan method of elimination is used, a test for zero divisor is 
included. A typewrlEeFmessage indicates when a pivotal element is smaller 

in absolute value than a value selected by the operator. At this point the solution 
may be continued or the next problem read in. 

C. Restrictions and Range : A maximum of 25 equations in 25 unknowns may be 
solved. A maximum of 8 significant digits per matrix element is allowed. 

D. Accuracy : Varies with the size of the matrix. 

E. Machine Configuration : Memory 20K, Card Input-Output. No other special 
features are required. 

F. Program Requirements : The entire core for 25 equations. 

G Source Language : The program is written in Fortran without Format. Alpha- 
betic output and format control are accomplished by two subroutines incorporated 
in the program. All computation is done in standard Fortran single precision 
floating point arithmetic. 

H. Program Execution Time: Read and compute time for 3 equations is 5 seconds. 
Typing the answer takes an additional 7 seconds. Read and compute time for 
8 equations is 25 seconds. 19 additional seconds are required for typing the 



I. Check Out Status : The program has been used successfully on systems of 3 to 
equations, and the logic has been verified for all other cases. 



answer. 




J. Sample Problem Running Time: N/A 

K. Comments: This program and its documentation were written by an IBM 
employee. It was developed for a specific purpose and submitted for generaL 
distribution to interested parties in the hope that it might prove helpful to 
other members of the data processing community. The program and its 
documentation are essentially in the author's original form. IBM serves as 
the distribution agency in supplying this program. Questions concerning the 
use of the program should be directed to the author's attention. 
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Description of the Program 



This program solves sets of non-homogeneous simultaneous linear (^) 
equations and provides either typewriter or punched card output. Up 
to 2.C equations in 2f unknowns may be solved. The first card read for 
each ¥eT^bT^quations state s the number of equations following. The 
program automatically sizes itself for the proper matrix: dimensions. 
Any number of sets of equations may be solved one after the other 
without reloading the program deck. When a solution has been obtained 
the word SOLUTION is typed. Under this heading the solution variables 
are typed one per line in the order in which they appear from left to 
right in the original equations. Punched card output is identical, with 
one card punched for each line typed. 



If a set of equations is submitted which has more than one solution, 477^ s 6iS ^ 



a valid solution will be typed or punched, but no indication will be 
given that there are other solutions as well. 



If the equations in the set are not all independent, a typewriter 
message will indicate that a diagonal e leme nt is ze r o. No solution is 
developed in this case, as there are infinite solutions to the set of 
equations. ~ " y ' v -z / ./>.,> ;> 

Additional information regarding the characteristics of the program 
is contained in the discussion of the method of computation. 



© 



Method of Computation 

General Description. The Jordan method of elimination* is used 
to accomplish the solution. Briefly this method consists of diagonalizing 
the augmented coefficient matrix. After diagonalization the right hand 
or augmented column contains the solution values for the variables. 

The reduction of the matrix takes place in the same storage locations 
where the original coefficients are stored. Thus at the completion of 
the solution the original coe fficients are los t, and in their place remain 
an identity matrix and the column of answers. 

Division by a Zero Diagonal . In the process of diagonalization the 
elements of each row to the right of the diagonal and the diagonal itself 
are divided by the diagonal element. The diagonal elements are used 
in the same order in which they appear in the original matrix. No attempt 
is made to select these pivotal elements by magnitude. Hence a zero 
must not appearjon Jihe_d^ If such a condition exists, 

a message will by typed indicating which diagonal element is zero and, 
control passes to reading the next problem. 

1 Alston S. Householder, Principles of Numerical Analysis (New York: 
McGraw-Hill Book Co. Inc., 1953), pp. 68-72. 



The possibility of zero diagonal elements may be minimized by 
observing two rules in laying out the original matrix. 

1. The upper left hand element should not be zero. 

2. No diagonal element may be zero if all elements in 
the row to the left of it are zero. 

There is a slight chance of generating a zero diagonal element in 
the course of calculations. Equations which are not independent (a 
unique solution undefined) will also generate a zero diagonal. If it is 
known that the coefficient matrix is non- singular, that is, that there 
is a unique solution, division by zero may Be~easily corrected and 
the problem submitted again. A simple rearrangement of the columns 
will usually suffice. In rearranging the columns "the e quality vector 
always must be kept on the right. As indicated later, rearrangement 
of columns may be accomplished simply by rearranging input cards. 
No additional punching is required. 

Tolerance Warning Level . When the diagonal element is very small 
relative to the values to its right in the row, a loss of accuracy will 
result from round- off errors. In this program the operator has control 
over this situation. When the absolute value of the divisor is less than or 
equal to some tolerance level selected by the operator, a warning message 
is typed. The operator then has the choice of continuing with the 
solution or proceeding to the next problem. If he should elect to 
continue, he does so with the knowledge that the answer may contain 
significant rounding errors. Of course, checking the answers will 
quickly show the accuracy of the solution. 



Input Card Format 



O 



The first card of the input for each set of equations must be a 
card with the number of equations punched anywhere between columns 1 
and 72 inclusive. If the number of equations is less than 10, the 
leading zero need not be present, although its presence makes no 
differenceT A decimal point must not appear with this number. The 
remainder of the card is blank. 

The matrix of coefficients is punched column-wise beginning with 
the left-most column and ending with the equality column. Each column 
must begin a new card. A card must contain five v alues. The column 
may be continued over as many cards as necessary. ffiffii& ca rds is the 
maximum number, as the program is limited to 2f) unknowns.) If 

the number of rows is not a multiple of five, dummy values must 

be punched to finish out five entries on the last card of the column. 
It is never necessary to punch more than four dummy entries for 

a given column. The number of cards per column will be the same for 
all columns within a given problem, but may vary from problem to 
problem. The dummy entries used to fill out the five entries on the 
last card of a column need not be zero. These entries are accessed 
during the read in, but are never accessed during calculation. Hence 

they must be present, but their value is immaterial. /f"^ 



The numerical value of each entry is punched in a free form notation- - 
the standard floating point variable form of 1620 FORTRAN. For those 
not familiar with this FORTRAN format, the specifics are detailed below. 

Each value may contain as many as 8 significant digits with a 
decimal point at the beginning, among, or~at the end of the group of 
digits. Zeros immediately to the right of the decimal point in a 
fractional number are not considered significant. A plus or minus 
sign may precede the number; the presence of the plus sign for , 
positive numbers is optional. Any value may be followed by an EC xx, 
where xx is a 1 or 2 digit power of 10 by which the number is to be 
multiplied. The plus (sigh^is optional, but if omitted the space it 
would have occupied must also be omitted. 

Examples of input values might be as follows: 



4-1. 



23.004678 -8.4 



.000000012345678 3.4E-A8 -.4E10 267.E-9 



In the last three examples the numbers represented are in order 
340000000. ,-4000000000. , and .000000267 . 



A DECIMAL POINT MUST APPEAR IN EACH ENTRY. 




At least one blank q plumn must be left between each of the five 
entries. The entries including leading blank columns may be in columns 
1-72 inclusive. 

Columns 73-80_pi each input card are not read by the program and 
are available for any identification the user may desire. 

Output Card Format 

The first card of output for each problem will be a header card with 
SOLUTION-^ punched in columns 1-10. Columns 77 -BO of this card will contain 
0413. The following cards contain the solution variables, one per card, 
in the order in which the variables appeared in the original matrix. The 
values are punched in standard 1620 FORTRAN floating point notation, 
that is, exactly like the input. However, the E - xx form will appear 
only if the decimal point does not fall at the beginning, end, or among 
the 8 significant digits. The punching begins in column 1. 

Output Typewriter Format 

This output is identical to the card output except that one line 
is printed for each card punched and the record mark following the 
word SOLUTION is not printed. This record mark stops transmission so 
the constant 0413 is not printed. 

Deck Identification 

The FORTRAN source language deck has SLE in columns 73-75 and 
a sequence number 01-59 in column 79-80. 

The machine language object deck in columns 77-80 contains a 
sequence number, u"00 1-0413. 



OPERATING INSTRUCTIONS 



I. A. Initial Console Setting 





Program 


Stop 


Parity 




X 


I/O 




X 


Overflow 


X 





B. Sense Switch Settings 



Sense Switch 


DN 


Solution is punched. 


1 


DFF 


Solution is typed. 


Sense Switch 
2 


DN 


Next problem after tolerance stop* 


DFF 


Continue after tolerance stop. 


Sense Switch 
3 


DN 


Tolerance level entered for each problem. 


DFF 


Tolerance level entered only once. 


Sense Switch 
4 


DN 


To correct error in typing tolerance level. 


DFF 


Tolerance level entered correctly. 



See Section V for further comment on Switches 2, 3, and 4. 



II. Input - Output 



Card Reader 



No. of Cards 


Description 


413 


Program Deck 


1 


N card - no of equations in 1st set of equations 


3-130 


Coefficient cards for 1st set of equations 


1 


N card - no of equations in 2nd set of equations 


3-130 


Coefficient cards for 2nd set of equations 




Etc. 
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Card Punch 



No. of Cards 


Description 


1 


SOLUTION Header card 


2-25 


Solution values for 1st set of equations 


1 


SOLUTION Header card 


Z-2£ 


Solution values for 2nd set of equations 




Etc. 


Typewriter Output 


Control 


Description 


Margins 


Set left margin as desired 


Tab Stops 


None 


Forms 


Standard paper 



III. Normal Loading Procedure 

1. Clear Storage 

2. Depress RESET 

3. Depress LOAD button 

4. Depress COMPUTER START after message LOAD DATA has been 
typed. See paragraph below. 

5. Depress READER START when hopper is empty 

When the LOAD DATA message appears the last few cards of the program 
deck will still be in the read hopper. If data cards are in the hopper 
following the program deck, pressing COMPUTER START will cause 
the remaining program cards to be read and a transfer to the program 
with subsequent reading of data cards. If no data follows the program, 
depressing COMPUTER START will cause loading to continue. The 
READER START must be depressed when the reader hopper is empty 
in order to read the last few program cards. When data cards are 
later placed in the hopper, READER START must be depressed to 
begin data entry and again when the hopper is empty to read the 
last data cards. 

IV. Special Loading Instructions 

At any time the program may be initialized and started by branching 
to 07548. 

6j 

V. Other Instructions and Remarks 

Tolerance Level. When the program begins, the typewriter will 
request the operator to enter the tolerance level as a floating point 
number, that is, in the same format as the input matrix coefficients- - 
see page 5. A record mark must immediately follow this entry. 



If the entry is typed correctly, set Sense Switch 4 OFF, press 

-8- 



RELEASE and COMPUTER START. The data cards will then be read. 



If the entry is typed incorrectly, turn Sense Switch 4 ON, press 
RELEASE and START. The program will immediately return control 
to the typewriter so the entry may be made again. If the retry 
is correct, follow the procedure for a correct entry. 

If it is desired to enter a tolerance level for each set of equations, 
turn Sense Switch 3 ON. If not, turn Switch 3 OFF, and the last tolerance 
level entered will automatically be applied to all succeeding sets 
of equations. 

When a tolerance warning stop occurs, directions for the use of 
Switch 2 are typed. Only at this time is the setting of Switch 2 
interrogated. If the tolerance level is greater than 9 integers or 
smaller than .000000001, the tolerance level will appear as an 
excess 50 floating point number in the message. 



All programmed stops are accompanied by typewriter messages 
which are self explanatory and which indicate the required action. 



As the program is written in FORTRAN, all the error stops 
and messages of that system apply. The more cryptic of these 
are of the form El, E2, E9, Fl, F2, F7. 

The FORTRAN processor without format, 1620-FO-002, was used 
to compile this program. All references to FORTRAN refer to this 
processor. 




VI. Programmed Stops and Required Action 





SAMPLE PROBLEMS 
The following four problems were used one after the other, 

8x1 ~ 5x 2 +7x 3 = 2 9 
+2x£ -2x3 = -2 
- Xj +9x2 -6x3 = 1 





Answer: x 


! = 2, 


x 2 = 


x l 


+2x 2 


+3x 3 


-4x 4 


+5x 5 




7x2 






-2x 5 


2x x 


+3x 2 


" 4x 3 


-5x 4 




X. 


+ X2 


+ X3 


+ x 4 


+ x 5 


lOxj 


+3x 2 


-4x3 




+ x 5 


3x x 


_3x 2 


-2x 


+2x 4 




-8x. x 




-9x 3 




+7x 5 



■ox^ r 1 X7 -oxg 

+ x ? + Xg 

+6x 6 

+ Xz + X*, + X 



+2x ( 

+ X/ 



-9x 7 

+6x7 
+3^ 



8 



-3x. 
-5x, 



Answers: xj = 1, = 2, X3 = 3, x 4 = 4 
x 5 = 5 » x 6 = 6 ' x 7 = 1 ' x 8 = 8 



-28 
19 
12 
36 

-42 
5 

18 
-19 



+2x- 



2x, +4x^ +7xo 
x^ +4x2 "*"^ x 3 



+10x 4 


+42x 5 


= -20 




+ x 4 


+ 5x c 


= 40 


, *■ • 


+ 4x 4 


+ 6x 5 


= -18 




+ 9x 5 


= 12 




- 2x 4 


+ 10x 5 


= 7 





This problem will cause diagonal 3 to be zero. It was included 
merely to illustrate the typewriter message. 



3xj +8x 2 +6 x 3 +10x 4 +42x 5 = -20 

+2x 2 + x 4 + 5x 5 = 40 

+ .OIX3 + 4x 4 + 6x 5 = -18 

2xj +4x 2 +7 X3 + 9x 5 = 12 

Xj +4x 2 +6 X3 - 2x 4 +10x 5 = 7 

This problem will make diagonal 3 go below a tolerance level of 
0. 1 . It was included merely to illustrate this typewriter message. 

The sample problem deck contains SP in columns 73-74; the problem 
number, 1, 2, 3, or 4, in column 77; and a sequence number in columns 
79-80. The sequence numbers are 01-05 for problem 1, 01-19 for problem 
2, and 01-07 for problems 3 and 4. 



~ /o - 



SAMPLE PROBLEMS 



INPUT 




3 



8. 0. -1. 0. .0. 

-5. 2. 9. 0. 0. 

7. -2. -6. 0. 0. 

29. -2. 1. 0. 0. 



8 

1. 0. 2. 1. 10. 
3. -8. 0. D. 0. 

2. 7. 3. 1. 3. 
-3. 0. 0. 0. 0. 

3. 0. -4. 1. -4. 
-2. -9. 0. 0. 0. 
-4. 0. -5. 1 • 0. 
2. 0. 0. 0. 0. 
5. -2. 0. 1. 1. 

0. 7. 0. 0. 0. 
-6. 0. 6. 1* 2. 

1. 0. 0. 0. 0. 
7. 1. 0. 1. -9. 
0. 6. 3. 0. 0. 

-a. i. J. l. o. 

o. -3. -5. o. 0. '4y 
-28, 19. 12. 36. -42. 

5. &18. -19. 0. 0. 



5 



3. 


0. 


0. 


2. 


1. 


8. 


2. 


0. 




4. 


6. 


. 


0. 


7. 


6. 



10. 1. 4. 0. -2. 
42. 5. 6.. 9. 10. 
-20. 40. -18. 12. 7. 



5 



3. 


0. 


0. 


2. 


1. 


8. 


2. 


0. 


4. 


4. 


6. 


0. 


l.E- 


•2 


7. 6; 


10. 


1. 


4. 


0. 


-2. 


42. 


5. 


"6. 


9 . 


10. 



-20. 40. -18. 12. 7. 




-7/ ' 



COMMENTS ON THE TYPEWRITER LOG 



FOR SAMPLE PROBLEMS 



On the typewriter log, which follows, several points should be noted. 
Sense Switch 3 was on at the beginning. Note that a tolerance level 
was entered for each of the first three problems. After the tolerance 
was entered for the third problem, Sense Switch 3 was turned off. 
When diagonal 3 went to zero and COMPUTER START was pressed to 
continue with the fourth problem, the solution began immediately without 
the entry of a tolerance level, as Switch 3 was off. Note the tolerance level 
(0. 1) in the warning message of the fourth problem is the last entry made-- 
that for problem 3. When this warning message appeared, Switch 2 was 
turned off and calculation continued to obtain a solution. 





SAMPLE PROBLEMS 
ft ■% 1 TYPEWRITER LOG 



LOAD DATA 



ENTER TOLERANCE LEVEL AS FP NO WITH RM 
THEN PRESS RELEASE AND START 
.0000000001* 



SOLUTION 
2.0000000 
3.0000000 
4. 0000000 



ENTER TOLERANCE LEVEL AS FP NO WITH RM 

THEN PRESS RELEASE AND START 

1.E-9* 



SOLUTION 
1 .0000065 
1 .9999976 
3.0000053 

3.9999963 
^. 9999916 
5.9999990 
7.0000022 
8.0000012 



ENTER TOLERANCE LEVEL AS FP NO WITH RM 

THEN PRESS RELEASE AND START 

.1* 



DIAGONAL 3 IS ZERO, PUSH START FOR NEXT PROBLEM- 



DIAGONAL 3 IS LESS THAN OR EQUAL TO .100000000 
TURN SWITCH 2 OFF TO CONTINUE, OR 
TURN SWITCH 2 ON FOR NEXT PROBLEM 
THEN PUSH START 



SOLUTION 
-6.7605150 

3*1.405115 
-7.0773320. _ 

5. 9^39480 
-6.9508359 



FLOW CHART 



Solution of n Simultaneous Linear Equations 
by the Jordan Method of Elimination 



8 




(tern* Ia> 
For M£S5 




Type M£ssf\&£> 

[To WTtrR 
TolERRA/cF 



REM) 




Dummy statement so that an 
alphabetic subroutine statement 
will not be the first executable 
statement 





1=1 



MESSAGES 



I 



Dirt** ftit 




-AT- 




11' 



c 






t , 














f 



Ty pe WARM 

.Fort Sw>TQH 2. 




[OFF] 



Switch 



oh 



C 
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!aJSE£T 



IS ZffcO 




c 
c 
c 



SOLUTION OF N SIMULTANEOUS LINEAR EQUATIONS 




10 D = SAY(2200.) 
= SAY f- -7201 .) 
ACCEPT , TOLER 

1 READ „ N 
N1 = M -f- % 
DO 2 J 1 , N1 
DO 2 I. - 1 , N, 5 

2 READ. AH, J), A(I+1,J), A(I+2,J), A(I+3,J), A(l+4,J) 
DO th f - 1 , N, 1 
DNO - / 
DUG - A( I , I ) 
IF ( D i AG ) 3, 20, k 

3 IF fPtAG + TOLER) 5, 19, 19 
k IF u ; u - TOLER) 19, 19, 5 

5 DO 6 J = I, Nl 

6 A( I , J) - A( I ,J) / D I AG 
K = 1 

9 IF (K~' ) 11 , 13, 1V 

11 FELM1 = A(K, I ) 
DO 12 J = I , N1 

12 A(K,J) . A(K,J) - FELMT * A( I , J) 

13 K = K + 1 
IF (K-N) 9, 9, ]k 

1 if CONTINUE 
J = N 1 

IF (SENSE SWITCH 1) 15, 17 

x 15 o « say(9104.) . T ; 

DO 16 1 - 1 , N 
16 PUNCH, A( I , J) 
GO T 21 
"/ 17 D -= SAY(2200.) 

D = SAY (-91 01 .) 
DO 18 I = 1, N 

18 PRINT, A( I ,J) 

21 IF (SENSE SWITCH 3) 10, 1 

19 D = SAY(2200.) 
D = SAY(-1 101 .) 
Y . 10001320. 

... Z * OUT (DNO) 

D . SAY(310V.) 
\ Y « 10001399a 
Z . OUT (TOLER) 
D = SAY (-4 301.) 
PAUSE 

_ IF (SENSE SWITCH 2) 21, 5 - — - 

\ 20 D » SAY(2200.) 
\ n « SAY (-11 01.) 
A Y = 10001320. 
Z = OUT (DNO) 
D * SAY(2101 .) 

... PAUSE 

GO TO 21 
END 



* Comments follow the program listing. . 



ALPHABETIC GONSTANT CARDS 
(Last 9 cards of program deck) 
DI AGONALZ 

15 ZERO* PUSH START FOR NEXT PROBLEMZ 

IS LESS THAN OR EQUAL TOZ 
TURN SWITCH 2 OFF TO CONTINUE* ORZ 
TURN SWITCH 2 ON FOR NEXT PROBLEMZ 
THEN PUSH STARTZ 

ENTER TOLERANCE LEVEL AS FP NO WITH RMZ 
THEN PRESS RELEASE AND STARTZ 
SOLUTION Z 



Z is a record mark, 0-2-8 punch. 



NOTES ON THE PROGRAM LISTING 



Two functions are used in the program. The function SAY is a 
subroutine for alphabetic output and is described in 1620 Library File No, 1. 6.006. 
The function OUT is a subroutine by R. F, Steinhart for format control 
of decimal points and is described in issue #38 of 1620 Technical 
Publications, Systems Engineering, IBM Midwestern Region, where it is 
referred to as the RFS subroutine. A description of this RFS subroutine 
is also included in 1620 Users Group Data, Distribution #1, July, 1961. 



The alphabetic constants for messages are the last 9 cards of the object 
program deck. They are not generated by the FORTRAN program, but 
are actually input data cards to the SAY function. They are read by 
the second FORTRAN statement, D = SAY (1950.), and therefore are 
placed immediately behind the object deck generated by FORTRAN, 
These 9 cards need be of no concern to the user, as they are read 
only once when the program is loaded. 

The FORTRAN processor, 1620-FO-002, used to compile this program was 
modified to include these two functions. 



END OF PROGRAM T0136 
START OF TABLE T2359 



UNSUBSCRIPTED VARIABLES 

T9999 EXP 

T9989 EXPF 

T9979 LOG 

T9969 LOGF 

T9959 SQR 

T99^9 SQRF 

T9939 SIN 

T99?-9 SINF. 

T9919 COS 

T9909 COSF 

T9899 ATN 

T9889 ATNF 

T9879 OUT 

T9869 OUTF 

T9859 ABS 

T9849 ABSF 

T9839 SAY 

T9829 SAYF 

T2799 K 

T2779 

T2719 TOLER 

T2699 N 

T2689 N1 

T2669 J 

T2659 I 

T2639 DNO 

T2629 D I AG 

T2529 FELMT 

T2429 Y 

T2409 Z 

SUBSCRIPTED VARIABLES 

T9819 T2809 
FLOATING POINT CONSTANTS 

T2769 1^19500000 

T2739 5422000000 

T2729 5472010000 

T?.489 1^91040000 

T2459 5491010000 

T2439 5411010000 

T2419 5810001320 

T2399 5431010000 

T2389 5810001399 

T2379 5443010000 

T2369 5421010000. 

FIXED POINT CONSTANTS 

T2789 0001 



END OF MAPPING 



Hji ) 



STATEMENT NUMBERS 

T2749 57548 0010 

12709 07632 0001 

12679 07704 5002 ^ 

T2649 09232 - 5014 _ - 

T2619 08324 0003 

T2609 09972 5020 - 

T2599 58392 5004 

T2589 58460 - 5005 

T2579 59700 5019 

T2569 58472 5006 - _ 

T2559 5868.8 5009 

T2549 58756 001 1 - .. .. .... 

T2539 59128 5013 

T2519 58852 5012 - 

T2509 59312 5015 

T2499 59476 - 001 7 - - 

T2479 59348 5016 

T2469 59680 5021 

T2449 59560 5018 

TEMPORARY ACCUMULATORS 

T2759 500 



Commeiits follow; 
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COMMENTS ON THE FORTRAN SYMBOL TABLE MAP 



The data in the symbol table map is defined as follows: 



END OF PROGRAM The first available location after che last program 

instruction- -excluding subroutines 

START OF TABLE The first available location preceding the expanded 

symbol table. 

UNSUBSCRIPTED VARIABLES 

Field address Variable name 

SUBSCRIPTED VARIABLES 

Field address Field address Variable name Number of 

first element last element rows 

FLOATING POINT CONSTANTS 

Field address Constant 
FIXED POINT CONSTANTS 

Field address Constant 
STATEMENT NUMBERS 

Field address Instruction address Statement 
of compiled of first instruction Number 

branch compiled for 

statement 

TEMPORARY ACCUMULATORS 

Field address Accumulator 

number 

Fields are addressed from their units position (high core position); 
instructions are addressed from their high order position (low core position). 



